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Abstract— In this study, the dispersion equation was formulated with the 
same procedure as the linear wave theory dispersion equation formulation 
procedure, but in this study, the total velocity potential was used. The 
dispersion equation obtained has the same form as the linear wave theory 
dispersion equation. There is a difference in its coefficients which causes 
the obtained wavelength to be half of the wavelength dispersion equation of 
the linear wave theory. 





I. INTRODUCTION 

The solution to the Laplace equation with the variable 
separation method is in the form of a superposition of two 
functions, namelycosine, andsine (Dean, 1991). In Dean 
(1991), the formulation of the linear wave theory 
dispersion equation is done by using a single component, 
namely the cosineonly. In this study, the velocity potential 
equation was carried out at the characteristic point where 
the value of the cosinefunction is the same as the value of 
the sine function. Therefore, the total velocity potential is 
in the form of cosine. The function obtained at this 
characteristic point will apply to thecosine andsine. 

The formulation of the dispersion equation was carried out 
with the same procedure as the formulation of the linear 
wave theory dispersion equation, which is to formulate the 
water level equation using the Bernoulli (Milne-Thomson, 
1960), carried out on the water surface. Then, with the 
water level equation obtained, the equation was formulated 
dispersion using linearized Kinematic Free Surface 
Boundary Condition. 


Considering that the same velocity potential was used and 
carried out with the same procedure, a dispersion equation 
with the same shape but a different coefficient was 
obtained. The dispersion equation was analyzed using the 
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wavenumber conservation equation. It was found that the 
dispersion equation obtained was not a function of water 
depth. 

In deep water, wave steepness was calculated using the 
maximum wave height from Wiegel (1949,1964). Wave 
steepness was compared with the critical wave steepness 
of Michell (1893) and Toffoli et all (2010). 

The dispersion equation in shallow water is formulated 
using the wave number and the energy conservation 
equation. With those equations, equations for the change in 
wave height and wavenumber due to water depth change 
were formulated. Even though it does not aim to develop a 
shoaling model, it produces a simple shoaling equation. 


Il. GENERAL SOLUTION OF LAPLACE’S 
EQUATION 
The Laplace’s equation(Milne-Thomson, 1960)on a two- 
dimensional axis system, namely(x,z)in the following 


form, 

ao ao 

ee ae 0....(1) 
p = ġ(x,z,t)is the velocity potential,xis the horizontal 
axis located on the still water surface;zis the vertical 
axis;tis time. The properties of the 


potential(Milne-Thomson (1960))are as follows, 


velocity 
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dp _ 
dx 
dp _ 
dz 
uis the velocity of the water particles in the -x horizontal 


direction. Meanwhile,wis the velocity of the particles in 
the -z vertical direction. 

The solution (1) is analytically carried out using a well- 
known method, namely the Variable Separation Method, 
where the solution is considered to be a multiplication of 
three functions(Bland, 1961), namely, 

d(x,z,t) = X(x) Z(z) T(t) ....(2) 


X(x)is a functionxonly, Z(z)is a functionzonly andT (t)is 
a functiontonly, which gives the general solution 

@(x, z,t) = (Acoskx + B sin kx) 

(Ce*? + De~**)sin(ot)....(3) 


This equation is a total general solution of Laplace's 
equation, which consists of two components, namely the 
cos kxand thesin kx components. 


2.1. Analysis of the Constants in the Solution of the 
Laplace Equation 

In (3), there are several constants, namely A, B, C, and D, 
which need to be formulated. These constants can be 
obtained by working on boundary conditions, where the 
Laplace equation is known as a boundary value problem 
where the specific solution was obtained by working on 
boundary conditions. 


Working on the kinematic bottom boundary condition on 
the flat bottom, as done by Dean (1991), was used to 
formulate the constants. The kinematic bottom boundary 
conditions on the flat bottom are: 

Wp = 0....(4) 

Using (3), the velocity of the water particles in the -z 
vertical direction is 


do 
Va = —k(Acoskx + B sin kx) 


(Ce — De-**)sin(ct) 
Then 


W_p, = —k(Acoskx + B sin kx) 
(Ce—* — De*")sin(ct) 


Substitution to (4), 

—k(Acoskx + B sin kx)(Ce7*" — De*")sin(ot) = 0 
This condition can be occured in 

Ce*h — Dekh = 0 

Where, 
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C=] Der onc (5) 

Substituting this equation into (3) will get the equation, 
@(x,z,t) = 2De*"(Acoskx + B sin kx) 

coshk(h + z) sin(ot) ...... (6) 


To get the constants A andB,it was done that on a wave as 
a whole there is only one single characteristic, whereA = 
B, 
@(x,z,t) = 2ADe*™"(coskx + sin kx) 
coshk(h + z) sin(ot) 
Defined G = 2ADe*", 
P(x, z,t) = G(coskx + sin kx) 
coshk(h + z) sin(ot)....(7) 


Proving that there is only one wave characteristic. Then, 
an analysis was carried out on G through the velocity 
equation, for example, the vertical velocity w. 
w(x, z,t) = —Gk(coskx + sin kx) 

sinh k(h + z) sin(ot) 


The particle velocitywhas units(m/sec)while the 
wavenumber khas units(m~t). Then, the unit ofG 
is(m. m/sec), indicating thatGis the rate of transfer of 
wave energy, which must be a single value, implying 
thatA = B. 

The formulation can be done with other procedures, 
namely (6) done at the characteristic point where cos kx = 
sin kx, then (6) becomes, 

@(x,z,t) = 2De*"(A + B) cos kx 


coshk(h + z)sinot...(8) 


For thecosine component. As for thesine component, 
@(x,z,t) = 2De*"(A + B) sin kx 
cosh k(h + z) sinot...(9) 


Both in (9) and (10), can be defined: 
G=2De A + B) ...... (10) 
Then, (8) and (9) become: 
(x, z,t) = G coskx cosh k(h + z) sin(ot) 
(11) 
@(x,z,t) = G sin kx cosh k(h + z) sin(ot) 
...(12) 
This proves the singularity of the value of G. However, it 
also indicates that the analysis of wave dynamics using 
(11) or (12)usesG which is a combination of two energies. 


2.2, The Wavenumber Conservation Equation 

In solving the Laplace equation with the variable 
separation method, there is no flat bottom assumption. The 
flat bottom assumption is only in the formulation of the 
constant. Furthermore, if (7) is done on a sloping bottom, 
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dk 


there will be values of = — , and 26 In the following 
dx dx dx 


: : dk 
section, an analysis of the an Was done. 
mf 


It has been mentioned that in solving the Laplace equation 
with the Variable Separation Method, it is assumed that the 
flow potential consists of 3 components, as presented in 
(2). 

In this equation, 

X(x) = coskx, isa function ofxonly, 

Z(z) = cosh k (h + z), is a function ofzonly 

T(t) = sinot, is a function oftonly. 


Given the nature of the function, it must be 





dZ(z) 
=0 
dx 
Given Z(z) = coshk (h + 2), then 
Wha ee 
dx 





ak(n+5) 7 
ad = 0... (14) 
Atz = 0, 
dkh 
ERSO ee (15) 


MI. FORMULATION OF THE DISPERSION 
EQUATION 
The dispersion equation will be formulated using velocity 
potential (11) and on a flat bottom. The velocity of the 
particle in the xhorizontal direction and the zvertical 
direction is, 
do 


== = Gksinkxcoshk(h + z)sinot ..(16)w = 
— = = —Gkcos kx sinh k(h + z) sin ot 
AT) 


1.1. Water Surface Equationn (x, t) 

To get the dispersion equation, first, the water surface 
equation is formulated (Dean (1991)). The equation was 
formulated using the Bernoulli equation which is carried 
out on the surface, namely 


— Sn 4 (u? + w2) + gn + Pa = C(#)..(18) 

The 7 index shows that the relevant variable is applied to 
the surface of the wave, g is the acceleration due to 
gravity, p is the pressure acting on the fluid particles, p is 
the density of the fluid, and C(t) is a constant that can be 
used zero for the periodic function (Dean (1991)). 


pnis the pressure on the surface, i.e., atmospheric pressure, 
by using atmospheric pressure as the reference 
pressure,p, = 0.After entering the dynamic free surface 


boundary conditionp, = 0, (18) divided byg 
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a 

gdt 2g 
The 2"! term is kinetic energy, while the 3" term is 
potential energy. At a small wave amplitude, the 2™4 term 


ue+weytn=O0 ee (19) 


will be much smaller than the 3" term. Therefore, it can be 
ignored, the Bernoulli equation becomes, 


d 
gdt 
The obtained water level equation is: 
do, 
E t = 
n(x, t) ae 


Substituting the flow potential equation, 


n(x,t) = eee cos kx cos a t...(20) 


Go coshk(h+n). 


For a periodic function, is a constant. Defined 


g 
p Gocoshk(h +n) 
g 
However, it has been mentioned that the use of a single 
velocity potential (11), implies that G is the sum of the 
energies of the two waves, therefore the wave amplitude 


equation must be divided by 2, 
_ Gocoshk(h+n) 


yom BIORE less (21) 


2g 
Ais the wave amplitude. But, G can not be calculated using 
this equation. G must be calculated using 

gA 
os 2ocoshk(h +n) 
This will be discussed in the next paper. 


The surface water level equation becomes 
n(x,t) = Acoskx cos o t....(22) 


3.2. Dispersion Equation 
The dispersion equation is formulated using the Kinematic 
Free Surface Boundary Condition (KFSBC), 


a, a 
Wy =e tu (23) 


: a 
For long waves and small amplitudes, upand- are small 


numbers. Thus, un Tbecomes a very small and negligible 
number, andKFSBC become, 


a 
Wy = G (24) 
Substituting (17) and (22)into (24) and the same terms on 
the left and right sides of the equation cancel each other 
out, the following equation is obtained 

20? = gk tanhk(h +n) ...... (25) 

In (25) there is anvariable whose value needs to be known. 
From the beginning, the formulation was carried out with a 
single flow potential, namely (11), where the formulation 
was carried out at the characteristic point wherecoskx = 
sin kx, on the time variable, also used the characteristic 
point wherecos ot = sin ø t. Thus, a characteristic point 
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that is complete is the point wherecoskx = sinkx == 
1 ; 1 
5 V2andcos ot=sinot= 5 v2. 


The solution obtained at this point is valued for 
bothsinandcos. At this point, the water level elevation of 


(22) is, 
cA 
= 


Then (25) becomes, 
20? = gktanhk (h + Ê) .....(26) 


Then, the wave number conservation equation (14) was 





reviewed. 
dtanhk (h + *) i 1 dk (h +4) E 
dx ~ cosh?k (h + $) dx 7 


It was found that 
A 
tanhk (h + 2) =c 27 


where c is a constantnumber. 


Equation (26) must apply to both deep and shallow water. 
In deep water 


A 
tanhk (h + 2) = 1.....(28) 
Then the constant c in (27) is 1. The dispersion equation 


becomes, 
207 =ghk....... (29) 


Therefore, the dispersion equation is only for deep water. 


Equation (29) gives a wavelength L = “half of the well- 
known linear wave theory wavelengths in the deepwater of 


the form 
of = gk... (30) 


In Table (1), the comparison between wavelengths of (29), 
L29, and wavelengths of (30), L39, whereLz9 = 0.5L30. 
The wave steepness of each wavelength is calculated using 
the maximum wave height ofWiegel (1949,1964), namely 


T2 
Ho max = zt. (31) 





15.62" ° 
Table.1: Comparison of wavelengths (29) and (30) 
T L30 L29 A max Ho max 
(sec) m) m) L30 L29 





6 56.207 | 28.104 | 0.026 0.052 
7 76.504 | 38.252 | 0.026 0.052 
8 99.924 | 49.962 | 0.026 0.052 
9 126.466 | 63.233 | 0.026 0.052 
10 156.131 | 78.066 | 0.026 0.052 
11 188.919 | 94.459 | 0.026 0.052 
12 224.829 | 112.414 | 0.026 0.052 
13 263.861 | 131.931 | 0.026 0.052 
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14 306.017 | 153.008 | 0.026 0.052 
15 351.295 | 175.647 | 0.026 0.052 


























Wave steepness (=)generated by the two dispersion 
equations in Table (1), is very small, much smaller than 


the critical wave steepness ofMichell’s criteria (1893), 
#0142 2c (32) 


-= 
And criteria for Toffoli et al (2010) 
#20170 an. (33) 


-= 
From this, although (29) produces a wavelength that is 
much smaller than (30), it still needs further development 
in order to achieve a wavelength that gives a critical wave 
steepness according to the criteria of Michell (1893)or 
Toffoli et al (2010). 


In order to obtain a critical wave steepness following the 
Mitchell or Toffoli criteria, the left side (29) is multiplied 
by a coefficienty, 

2yo? = gk....(34) 


yis a coefficient that is greater than 1. By trial and error, 
the value of y = 2.75is obtained which produces a 
wavelength with a wave steepness corresponding to the 
critical wave steepness of Michell (1893), with the 
calculation results presented in Table (2). 
Table.2: Wavelength and wave steepness at 
y = 2.75 
T L Ho,max 
(sec) (m) | I 
6| 10.219 0.142 
7 13.91 0.142 
8 18.168 0.142 
9 | 22.994 0.142 
10 | 28.387 0.142 
11 | 34.349 0.142 
12 | 40.878 0.142 
13 | 47.975 0.142 
14 | 55.639 0.142 
15 | 63.872 0.142 





















































Meanwhile, if Toffoli et al criteria are used, it will 
obtainy = 3.285. 


However, (29) and (34), need more intensive analytical 
study , to fulfill the kinematic free surface boundary 


condition. 


3.3. Dispersion equation in shallow water 
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To obtain the dispersion equation in shallow water, the 
wave number conservation equation is carried out. First, 


define the value of kọ (ho +4) in deep water where 


tanhkko (ho + £) =1, 


ko (ho +€) = Or.......35) 


where tanh(@7) = 1. Index 0 indicates the variable in 
deep water. The value of@is a positive number whose 
value is greater than or equal to 1. SPM (1984) uses the 
value of = 1. According to the law of conservation of 
wavenumbers (14), 


An Ap 
kn(h+ >) 2 ko (ho +) = 6n 


Then 





The index h indicates the variable at the water depth h, 
smaller than hg. In this equation, there are two unknowns, 
namely kp and Ap. One more equation of the relation 
between kp and Ap is needed. The available equation is the 
energy conservation equation. The wave energy at one 
wavelength (Dean (1991)) is as follows, 


E =~pg H*L......37) 

g is the gravitational force, and p is the mass density of 
water. Assuming there is no loss of wave energy, then the 
relationship should be as follows, 

Holy, = HêLo.....(38) 

where Hp and Lo are wave height and wavelength in deep 
water, while H, and L, are wave height and wavelength at 
shallower water depth h. Equation (38) can be expressed 
as the equation for the wave number kp, 

ky = ath (39) 

Where Ey = HiLo, substitution (39) to (36) assuming a 
sinusoidal wave where A; = 0.5H,, the equation for H, is 
obtained. 


3 
fi + nH? — “2 = 0....(40) 


H, is calculated by (40), then kp is calculated by (39). 

For example, the calculation of wavelength in shallow 
water used waves with a wave period of 8 sec. deep water 
wave height Hy is calculated by (31), deep water 
wavelength Ly is calculated by (29) or by (34) with y = 1. 
The calculation is done by using 0 = 1, and ignoring 
breaking, the calculation results are compared with the 
linear wave theory dispersion equation, namely, 
o?=gktanhkh a... (41) 

The calculation results are presented in Fig (1) below. 
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100 
80 
GB 60 
— 40 
20 7 {ter | 
0 
0 5 10 15 20 25 


water depth h (m) 





eq. (41) 





eq. (39) 


Fig.l: Wavelength of (41) and (39). 


In Fig (2), the change in wave height to water depth 
(shoaling) is presented in the form of a nonlinear line. In 


the figure, the calculation is stopped when Z = 0.80 which 


is considered breaking. At that point, the wave height is 
H = 4.80 m, the water depth is h= 5.99 m, and the 
wavelength is L = 14.40 m. From the breaking conditions, 
it was found that the shoaling that occurred was too large. 


5 
4 


3 


Wave height H (m) 


0 5 10 15 20 25 
Water depth h (m) 


Fig.2: Shoaling on a wave period of § sec. 


The magnitude of the resulting shoaling is because the 
wave energy is too large. After all, the wavelength Ly 
calculated by (29) is too long. When used (34) using y = 
3.285 where the value of y is the result of critical wave 
steepness adjustment with the criteria of Toffoli et al 
(2010), it is obtained that the deep-water wavelength Lo is 
shorter. Thus, the wave energy decrease, and the resulting 
shoaling is also not too large. 


In Table (3), the results of the shoaling calculation are 
presented where Lọ is calculated by (34) using y = 
3.285,and water depth coefficient 0 = 1.The wave height 
at T= 0.78 was compared with the average breaker 
height of the five-breaker height index (BHI) from a 
number of previous studies conducted by Komar and 
Gaughan (1972), Larson, M. and Kraus, N.C. (1989), 
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Smith and Kraus (1990), Gourlay (1992), and Rattana 
Pitikon and Shibayama (2000). 
Komar and Gaughan (1972) 

1 


Hb _ Ho TE 

2 = 0.56 (=) re (42) 

Larson and Kraus (1989), 

Hp Ha) 70-24 

2 = 0.53 (=) ees (43) 

Smith and Kraus (1990), 

Hp Hy\70:30+0.88m 

i = (0.34 + 2.74m) (=) seins (44) 

Gourlay (1992), 

Hp _ H) 70-28 

72 = 0.478 (=) ene (45) 

Rattana Pitikon and Shibayama (2000) : 

1 

Hy Ho\ 5 
— = (10.02m? — 7.46m? + 1.32m + 0.55) (=) 
Ho Lo 

TEETE (46) 


In these BHI equations, Hois deep water wave height, Lo is 
deep water wavelength (calculated using linear wave 


2 2 ; : 
theory, ky = ToLo = 5, m is the bottom slope and H, is 
0 


breaker height. In this study, bottom slope m = 0 is used. 


Breaker depth h, in Table (3) is calculated by the breaker 
depth equation from SPM (1984), that is, 


hp 1 Hp 
aE 4H rhp = — p (47) 
Ho b-(3r2) p-(gr2) 
1.56 
a= 43.75(1 = e 190m) 5 = T+ e 195m 
Calculation with bottom slope m = 0 obtained Ze = 0.78. 
b 


Table.3: Comparison of equations (39+40) with 

















BHI. 
Eq. (39+40) BHI 
T Hn h H, hp 
(sec) (m) (m) (m) (m) 
6 1.805 | 2.314 | 1.721 | 2.207 
7 2.457 | 3.149 | 2.343 | 3.003 
8 3.209 | 4.114 3.06 | 3.923 





9 4.061 5.206 3.873 4.965 
10 5.014 6.427 4.781 6.129 
11 6.067 7.777 5.785 7.417 
12 1.22 9.256 6.885 8.826 
13 8.473 10.863 8.08 10.359 
14 9.827 12.598 9.371 12.014 
15 11.281 | 14.462 10.757 | 13.791 
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Therefore, Table (3) mentioned the value of = 2 = 
b 
0.78. It is found that H is close to H, and h is also close 
to hp. This indicates that the shoaling of (39+40) can be 
improved by shortening the deep wáter wavelength and the 


wavelength Ly calculated by (29) is still too long. 


IV. CONCLUSION 
A review of the critical wave steepness in deep water 
results that the wavelengths obtained using the total 
velocity potential are better than the wavelengths 
formulated using a single velocity potential, likewise with 
the results of studies in shallow water. However, it still 
needs to be shortened again. 
In short, it can be done in a simple way, by increasing the 
coefficients in the dispersion equation. However, it is 
necessary to investigate the origin of these coefficients 
analytically based on hydrodynamic equations, especially 
the kinematic free surface boundary condition and the 
Bernoulli equation. 
Considering that the formulation of the dispersion equation 
in deep water in this study was carried out following the 
procedure for formulating a linear wave dispersion 
equation and producing a dispersion equation in the same 
form as the linear wave theory dispersion equation, the 
dispersion equation obtained can also be referred to as a 
linear wave dispersion equation. 
To conclude, the analysis of wave dynamics using the 
velocity potential solution of the Laplace equation with the 
Variable Separation Method should use the total velocity 
potential. 
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